Evidence for existence of many pure ground states in 3d ±J Spin Glasses 
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Ground states of 3d EA Ising spin glasses are calculated for 
sizes up to 14 3 using a combination of genetic algorithms and 
cluster-exact approximation . The distribution P(\q\) of over- 
laps is calculated. For increasing size the width of P(|g|) con- 
verges to a nonzero value, indicating that many pure ground 
states exist for short range Ising spin glasses. 
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Introduction The behavior of the Edwards-Anderson 
(EA) ± J Ising spin glass with short range (i.e. realistic) 
interactions is still not well understood. The EA Ising 
spin glass is a system of N spins er { ; = ±1, described by 
the Hamiltonian 



H 



(i,3) 



0~i<Jj 



(1) 



The sum goes over nearest neighbors. In this let- 

ter we consider 3d cubic systems with periodic boundary 
conditions, N — I? spins and the exchange interactions 
(bonds) take = ±1 with equal probability under the 
constraint J2<ij> J*j = 0- 

One of the most important questions is, whether many 
pure states exist for realistic spin glasses. For the infi- 
nite ranged Sherrington-Kirkpatrik (SK) Ising spin glass 
this question was answered positively by the replica- 
symmetry-breaking mean- field (MF) scheme by Parisi |Q . 
Some authors g] disbelieve that the MF theory is valid 
in total for realistic spin glasses. But also a complete 
different model is proposed: the Droplet Scaling (DS) 
theory suggests that only two pure states (related 
by a global flip) exist and that the most relevant exci- 
tations are obtained by reversing large domains of spins 
(the droplets). In || ground states were calculated using 
multicanonical sampling, but no decision between MF 
and DS could be made. Recent results of Monte Carlo 
(MC) simulations || at finite temperature find evidence 
for the MF picture. 

In this letter we show our results for the calculation of 
spin glass ground states using a hybrid of genetic algo- 
rithms jH],[ll| and cluster-exact approximation (CEA) 
[ p2[ . CEA is a discrete optimization procedure which is 
based on the construction of spin clusters exhibiting no 
frustration and calculating exact ground states for these 



clusters using graph theoretical methods. By calculating 
the distribution of overlaps for different sizes L we find 
evidence for the existence of many pure ground states for 
the ± J spin glass. 

Observables For a fixed realization J — {Jij} of the ex- 
change interactions and two statistical independent con- 
figurations (replicas) {erf}, {erf}, the overlap |l[] is de- 
fined as 



a/3 — 



(2) 



The ground state of a given realization is characterized 
by the probability density Pj(q). Averaging over the 
realizations J, denoted by [■]„„, results in (Z = number 
of realizations) 



P iq ) = [ Pj ( q )} av = ±-Y, p j(<i) 



(3) 



Because no external field is present the densities are sym- 
metric: Pj(q) = Pj(—q) and P(q) = P(—q). So only 
averages of \q\ n are relevant: 



\qj I 



\q\ n Pj{q)dq 



-1 

1 



\q\ n P(q)dq 



(4) 
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The Droplet model predicts that only two pure states 
exist, implying that P(q) converges for L — > 00 to 
P(q) = \{8{q — \q\) + 5(q + \q\)) (we don't indicate the 
L dependence by an index), while in the MF picture the 
density remains nonzero for a range —q\<q<q\ with 
peaks at ±c/ max (0 < q max < <Zi)- Consequently the vari- 
ance 



\\q\) 



(M - \q\) 2 P(q)dq=\q\ 2 - \q\ 



(6) 



stays finite for L — > 00 in the MF pictures (ct 2 (|<7| 



o 2 / L 3 ) while ct 2 



according the DS 



ai 

approach. Here y is the zero-temperature scaling expo- 
nent @ denoted 6 in |,|. 

We have measured the energy per spin e as well. Addi- 
tionally to the usual fits other scaling corrections were 
discussed [^0[, including an exponential one (e = + 
ai exp(-a 2 L)). 
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Results We performed ground state calculations for 
sizes L = 3,4,5,6,8,10,12,14. We calculated from 100 
realizations for L = 14 up to 4000 realizations for L = 3. 
One L — 14 run needs typically 540 CPU-min on a 
80MHz PPC601 processor (70 CPU-min for L = 12, 
0.2 CPU-sec for L = 3). For each realization 40 indepen- 
dent runs were made. Each run resulted in one configu- 
ration, which was stored, if it exhibited the ground state 
energy. For L = 14 this resulted in an average of 9.9 
ground states per realization while for L = 3 we got an 
average of 39.9 ground states. Details of the algorithm, 
simulation parameters and more results are given in [|l3| . 
The results of the calculations are summarized in table 

Because we want to estimate the quality of the mini- 
mization procedure, we discuss the average ground state 
energy at first. By comparison with previous calcula- 
tions ||[l0) we can see that we have obtained the same 
or slightly lower (i.e. better) energy, except for L = 12 in 
|^| which is based only on 7 realizations. So we are con- 
fident that we get very close to the true ground states 
using our approximation procedure. 
We extrapolate to L — > oo by using finite size scaling 
(FSS) fits. We used the Levenberg-Marquard Method 
implemented in |l4|] which gives also the probability Q of 
finding such a bad data-set (provided that the fit function 
is exact). 

Table [n] show the results for FFS fits for the following 
three functions: 

f(L) = + ai/L 3 (7) 
g(L) = aoo + ai L- a * (8) 
h{L) = aoo + «i exp(— a2-L) (9) 

For cr 2 (|q|) results for fits of the first two functions f{L) 
and g(L) (with a^ set to 0) are given in the same table 
as well. 

We get a ground state energy of the infinite system of 
Goo = aoo = -1.7868(3) (Q=0.20) resp. 6^ = 0^ = 
-1.7876(3) with a 2 = 2.84(3) (Q=0.77). The exponential 
fit is ruled out (Q=10~ 7 ). Our ground state energy is 
lower than some previous results: = —1.75 |HJ, = 
— 1.76(2) and of the same magnitude as the most 
accurate study so far = —1.7863(4) flPfl . In || a 
lower value was found (e^ = —1.839(4)), but that study 
suffers from only seven L = 12 realizations. 
The result for the energy, together with the best fit, is 
displayed in fig. |l|. 

In [[l0| a L~ 3 fit was ruled out, instead a L~ 2 97 fit was 
preferred. Our data shows a similar tendency, but the 
data for the energy is not accurate enough to decide this 
question. The second fit supports the MF picture too, 
because the value y = a 2 = 2.84 of the exponent exceeds 
the limit y < (d— l)/2 proposed by DS. Since the aim of 
this letter is to distinguish whether only two or infinitely 



many pure ground states exist, we proceed by analyzing 
the order parameter distribution P(\q\). 
In fig. H the probability density of the overlap is displayed 
for L = 6 and L = 12. Each realization enters the sum 
with the same weight, even if different numbers of ground 
states were available for the calculation of Pj(\q\). Only 
a small difference for high q-values is visible, which is 
reflected also in a slight shift of the average to smaller 
values (see table ^). But especially for \q\ < 0.8 no sig- 
nificant reduction in the probability is visible. So many 
different ground states still exist for higher L. This re- 
sult is supported quantitatively by the variances of \q\ as 
function of L which are shown in fig. |^. For all sizes the 
magnitude of the values is higher than in indicating 
that their MC algorithm does not explore the configura- 
tion space sufficiently. Our fits show a higher probability 
for the existence of many pure states (Q — 0.79) against 
the two state picture (Q = 0.09). 

Another interesting questions is whether the difference 
between the density functions of single realizations and 
the average -P(|g|) remains nonzero for increasing L, i.e. 
whether the short range spin glasses exhibit lack of self- 
averaging. But presently we have not computed enough 
ground states per realization for large sizes to address 
this question. 

Conclusion By the calculation of ground states using 
genetic cluster-exact approximation we find evidence for 
the existence of many pure states in short range spin 
glasses. The major drawback of our method is that the 
ground states are not exact, but for these lattice sizes this 
is beyond the power of actual algorithms and computers. 
By comparison with other results we feel confident that 
the energy of our states is very close to the true ground 
state energy. In the moment we do not have enough 
ground states per realization calculated (for higher L) to 
decide whether short range spin glasses exhibit lack of 
self-averaging or not. 
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FIG. 2. P(\q\) for £, = 6,12 
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FIG. 1. Ground state energy for 3 < L < 14 and best FSS 



fit 
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FIG. 3. Variances a 2 (\q\) for 3 < L < 14 and FSS fits 



3 



L e \q\ a 2 (\g\) 

3 -1.6731(19) 0.742( 3) 0.0731( 9) 

4 -1.7370( 9) 0.748( 3) 0.0669( 9) 

5 -1.7603( 8) 0.743( 4) 0.0628(13) 

6 -1.7723( 7) 0.742( 4) 0.0632(15) 
8 -1.7802( 5) 0.726( 6) 0.0621(20) 
10 -1.7840( 4) 0.737( 6) 0.0588(22) 
12 -1.7851( 4) 0.698( 8) 0.0593(30) 
14 -1.7865( 7) 0.701(20) 0.0583(72) 



TABLE I. Results. 



fit ai ci2 Q 

e f(L) -1.7869(3) 3.15( 5) - O20 

e g(L) -1.7876(3) 2.60(15) 2.84(5) 0.77 

e h(L) -1.7846(8) 1.01(17) 0.75(5) 10~ 7 

a 2 (|g|) f(L) 0.0608(6) 0.34(3) - 0.79 

a 2 (\q\) g(L) 0.088(4) 0.18(3) 0.09 



TABLE II. Result of FSS-fits. 
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